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Cn ' Abstract. We investigate the three-dimensional structure of the nearby edge-on spiral galaxy 
NGC 891 using 3D Monte Carlo radiative transfer models, with realistic spiral structure and 
fractally clumped dust. Using the spiral and dumpiness parameters found from recently com- 
pleted scattered hght models we produce lower resolution SED models which reproduce the 
global UV-to-FIR SED of NGC 891. Our models contain a color gradient across the major axis 
(-H of the galaxy - similar to what is seen in images of the NGC 891. With minor adjustment our 

O ,' SED models are able to match the majority of M51s SED, a similar galaxy at a near face-on 
I different inclination. 
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1. Introduction 

^ • The interplay between starlight and dust in spiral galaxies is enormously complex. In 

the last decade, space-based observatories have allowed astronomers to measure spectral 



S^ ■ energy distributions (SEDs) from the ultraviolet (UV) through the far-infrared (FIR) 
f^ I These SEDs show that radiative transfer (RT) models based only on optical data (e.g. 
^^ ' Xilouris et al. 1999) significantly overpredict the amount of escaping UV radiation and 
^— s ^ underpredict the amount of dust present in edge-on spirals (Popescu et al. 2000; Bianchi 
p^ ■ 2008, hereafter BOS). To combat this problem astronomers have created more complex 
y—i ' models, including additional stellar disk components and semi-analytic approximations 
^ . for the effects of enshrouded star formation with clumpy dust and light. 
• rH I Most of these advanced models have focused exclusively on fitting galaxy SEDs at the 

K% ■ expense of creating models with observationally accurate images, due to the difficulty in 
S ' fitting non-axisymmetric parameters like clumpy dust substructures or spiral arms. 
" " ■ We have recently completed a preliminary investigation into optimizing non-axisymmetric 

RT models fit to high-resolution Hubble Space Telescope optical images of the central re- 
gion of the edge-on spiral galaxy NGC 891 (Schechtman-Rook et al. 2011). These models 
included fractally clumped dust and logarithmic spiral arms specially designed to match 
observations of the arm/interarm ratios of face-on spiral galaxies. The best-fitting mod- 
els show a strong preference for spirality and dumpiness. Here we investigate the effect 
of the spirality and clumping parameters determined from these fits to optical images 
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Table 1. Model Parameters 
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on fitting tlic SEDs of NGC 891 and M51. We use a new 3D Monte Carlo RT code 
(HYPERION, Robitaille et al. 2011), which uses arbitrary dust and emissivity grids and 
includes multiple dust grain sizes (including PAHs). 

2. Model 

Our 100x100x50 pixel SED models are mapped to a 40x40x20 kpc volume. Stellar 
emission is set to the Sb spectrum of Fioc & Rocca-Volmerange (1997), split into a old 
population of stars with masses <3 M0(bulge and Disk-1) and an young population 
(M > 3 Mq; Disk-2). Model parameters are shown in Table 1. Md is the total dust 
mass, Li and L2 are the total bolometric disk luminosities, and T denotes the density 
threshold for young stellar emission (where appropriate). All models have a Sersic n=4 
bulge with Lbuigc — 1-8 x lO^^L© and structural parameters identical to those used in 
Bianchi (2008). The model ratios of old-to-young stellar populations span unredenned, 
spatially-integrated SEDs corresponding to a wide range of Hubble types (for reference 
template E,Sa,Sc,Im galaxies from Fioc & Rocca-Volmerange (1997) have old-to-young 
ratios of 35,9.8,2.3,1.1, respectively). 

Disk-1 has half the scale-length and twice the scale- height of the dust. Disk-2 has the 
same scales as the dust, and when it is clumpy it follows the distribution of the clumpy 
dust; higher dust density cells weighted to have more stellar emission, dumpiness and 
spirality are consistent with what we found in Schechtman-Rook et al. (2011), i.e., 50% 
of the dust mass is put in clumps distributed in 130 parent cells; spiral arms have a pitch 
angle of 20° and a relative strength of 50%. For models 4 and 5 there is a minimum dust 
density required to place Disk-2 emission in a cell; this threshold mimics the effects of 
embedded star formation. 



3. Results 

Single disk models are unable to simultaneously match the UV and FIR SED of NGC 
891, even when the dust mass and stellar luminosity are increased or dumpiness and 
spirality are added. The poor fit from single disk models was also noted by B08. Adding 
spirality, regardless of the amplitude of the spiral perturbations, has almost no effect on 
the spatially-integrated SED morphology, a result found also by Popescu et al. (2011). 

Shown in Figure [l] restricting the young stellar emission to areas of higher dust density 
(Model 4) allows us to dramatically increase the luminosity of the second disk without 
increasing the dust mass. This density threshold appears to be unique for any combination 
of disk luminosities and dust masses — higher thresholds push the dust peak too far into 
the blue, while at lower thresholds the model cannot simultaneously reproduce the UV, 
FIR, and 12 /im PAH feature. While this model still predicts a low 24 /xm fiux, high- 
resolution simulations of individual dust clouds indicate that this discrepancy is due to 
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Figure 1. Spatially integrated SEDs. Left: Two-disk models (lines) compared to data from 
NGC 891 (red points, from BOS). Pink and blue shading denote contributions to Model 4's SED 
from old and young stars, respectively. Right: Face-on two-disk models compared to data for 
M51 (from the NASA/IPAC ExtragalacUc Database). 

the unresolved nature of the dense clumps in our simulations (BOS). We investigated 
the viability of our models at other projected inclinations by comparing them to M51 
(i ~ 20deg; Figure 1). For simplicity we fix all geometric parameters (including clumping 
and spirality) of the dust and stellar distributions as well as L2, and allow only Mj and 
Li to vary. A good match to the observed, spatially integrated SED (except at 24/im, 
as explained above) is found by increasing the dust mass by 75% and halving the old- 
disk luminosity (model 5). The FIR flux distributions indicate that M51 has a dearth of 
cold dust compared to N891 despite M51's increased dust mass, consistent with elevated 
star- formation induced by M51's ongoing interaction. 

While the addition of spirality has little effect on the shape of the spatially-integrated 
SED, it does produce a gradient in the UV/blue light across the major axis of the 
model. This gradient is present in NGC 891 and is commonly considered to be due to 



the presence of a a grand design spiral pattern Kamphuis et al. (2007) This gradient 



presents an opportunity to better constrain the otherwise hidden thin disk of edge-ons, 
and can only be fit by models with non-axisymmetric structure like ours. 
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